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We report electronic transport of Dirac cones when Fe is replaced by Ru, which has an isoelectronic 
electron configuration to Fe, using single crystals of Ba(Fei_ :r Ru :E As)2. The electronic transport 
of parabolic bands is shown to be suppressed by scattering due to the crystal lattice distortion 
and the impurity effect of Ru, while that of the Dirac cone is not significantly reduced due to the 
intrinsic character of Dirac cones. It is clearly shown from magnetoresistance and Hall coefficient 
measurements that the inverse of average mobility, proportional to cyclotron effective mass, develops 
as the square root of the carrier number (n) of the Dirac cones. This is the unique character of the 
Dirac cone linear dispersion relationship. Scattering of Ru on the Dirac cones is discussed in terms 
of the estimated mean free path using experimental parameters. 

PACS numbers: 74.70.Xa, 74.25.Dw, 72.15.Gd, 75.47.-m 



I. INTRODUCTION 

Transport phenomena of massless fermions are intrigu- 
ing subjects in the physics of Dirac cones and have be- 
come a very important frontier research area. A Dirac 
cone characterized by a degeneracy point in the linear 
band dispersion is known to be realized in materials with 
special geometrical symmetries, such as the K-K' degen- 
erate points in graphene [l|, topological insulators (TIs) 
with the Z? symmetrical constraint |2( , some organic con- 
ductors Q with a degenerate crossing band, and iron 
pnictides with different d^-d^ orbital symmetry (4T4l6l|. 
One of the significant features of electronic transport in 
the Dirac cone states is the large suppression of backward 
scattering as a consequence of the Berry's phase. 

In iron pnictides, the Dirac cone is created in pairs as 
the node of the spin density wave (SDW) ordering states 
caused by the orbital symmetry in the Fe-multi band sys- 
tem p| . The resulting Dirac cones have pseudo-spin chi- 
rality caused by the chirality compensation between the 
Dirac cones and the parabolic hole pockets @. These 
Dirac cones might therefore cause a large suppression in 
backward scattering. The situation can be compared to 
that observed for graphene, where the unit cell has two 
sublattices generating the Dirac cones in pairs at points 
K and K' in the first Brillouin Zone with opposite chi- 
rarities p|. The absence of backward scattering holds 
for the geometrically identical K and K' points, but it 
is broken in the case of intervalley scattering due to the 
opposite chirarity [13] ■ In recent studies on topological 
insulators, a single Dirac cone emerges to create a heli- 
cal spin texture at the time reversal invariant momentum 
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due to spin orbital coupling, and therefore the locked spin 
angular momentum prohibits backward scattering Q. 

In iron pnictides, almost no superconducting fluctua- 
tion has been reported in the underdoped regime of iron 
pnictide superconductors [l8j . which is in strong contrast 
with what has been encountered for the high tempera- 
ture superconducting cuprates flfj| - [23| . If the high mo- 
bility carrier could exist behind the Cooper paired elec- 
tronic states, the phase coherence of Cooper pairs would 
be much enhanced, which can result in a robust super- 
conducting state against thermal fluctuations. This is 
significantly important in the physics of the high tem- 
perature superconducting materials. 

Recently, a pronounced contribution of parabolic 
bands to electronic transport has been reported in 
Shubnikov-de Haas oscillation and Hall coefficient mea- 
surements for high quality single crystals made by an- 
nealing under BaAsj2J, [25|. Recent magnetoresistance 
measurements have reported that the Dirac cone per- 
sists in the underdoped regime of Ba(Fei_ a; Ru a ;As)2 [2a |. 
Considering the present situation established in the Dirac 
cone states in Ba(Fei_ 2; Ru 2 ;As)2, it is very important to 
clarify the scattering effect of the Dirac cone states of this 
family from the viewpoint of electric transport in detail, 
although the contribution of the parabolic bands must 
be suppressed to extract the intrinsically unique carrier 
transport of the Dirac cone in a complex multi band sys- 
tem of iron pnictide materials. In this research, we have 
studied this intriguing electronic state in the Fe-multi 
band system using Ru-substituted Ba(Fei_ x Ru x As)2- 
For this purpose, single crystals of Ba(Fei_ a: Ru :c As)2 
were rapidly quenched at around 1050 - 1100°C, which 
can introduce a significant lattice distortion effect as dis- 
cussed in detail in the text. It will be shown that the 
mobility is greatly reduced by scattering due to crystal 
lattice distortion in the case of the parabolic bands, while 
that of the Dirac cone is not significantly influenced as 
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the loss of backward scattering. We will show that the 
inverse of the average mobility (l/// a vc)j proportional to 
the cyclotron effective mass (m*), evolves as a function 
of square root of the carrier number, n, of the Dirac cone. 
The n dependence of m* is the unique character of the 
massless Dirac quasiparticles. It will be shown that the 
estimated mean free path is constant for the whole region 
of x and is not correlated with the Ru-Ru distance. 



II. EXPERIMENTAL 



Single crystals of Ba(Fei r Ru r As)2 were grown by a 
flux method using FeAs [2a, [27| . The quality of the sin- 
gle crystals was studied by synchrotron X-ray diffraction 
measurements at the beam line BL02B2, SPring-8. Elec- 
tric resistivity p measurements were also carried out using 
a four-probe method to check the quality of the samples. 
The Ru concentration x for Ba(Fei_ a; Ru 2 ;As)2 was deter- 
mined by employing the relationship between the c-axis 
lattice constant and x as reported in Ba(Fei_ 2: Ru l; As)2 
[26l [2?| ■ The dependence of in-plane transverse magne- 
toresistance (Rm) and the Hall coefficient (Ra) on mag- 
netic field B was measured using a four-probe method in 
B of -9 - 9 T at various temperatures (T) between 2 and 
300 K. 



S 2 

u 

a 



1 t ' 

w 


1 

• x = 




□ 0.022 




♦ 0.051 

0.102 






T = 2K 


— 



i 



(b) & 


I I I I I I 


r = 2K 


f nOo° 000000 ° 
► o 0°° 


00000000 00 «J 
♦♦*T 


□ 0.022 
♦ 0.051 
0.102 
■ 0.150 



-8 -4 4 8 

B(T) 

FIG. 1: (color online) Magnetic field (B) dependence of the 
Hall resistivity (p yx ) (a) and the derivative of magnetoresis- 
tance (di?M/dB) (b) for Ba(Fei_ a; Ru a ;As)2 with x between 
and 0.150 at 2 K. 



III. RESULTS AND DISCUSSIONS 



as 



Figures 1 (a) and (b) show the B dependencies of 
Hall resistivity (p yx ) and the derivatives of magnctoresis- 
tances (dRu/dB) for Ba(Fei_ 2; Ru 2 ;As)2 with x between 
and 0.150 at 2 K. For x = 0, p yx developed linearly with 
B in the low B regime, but showed deviation from the an- 
ticipated straight line. This is consistent with the nonlin- 
ear behavior of p yx in the multiband nature of Ba(FcAs)2 
p8l [29| . The nonlinear behavior of p yx with B was sup- 
pressed by Ru doping and an almost linear p yx was ob- 
served at x = 0.150. The value of dRyi/dB developed 
linearly with B but was saturated above 2 T for x = 0, 
which indicates that the B dependence of Rm from Rm 
oc B 2 to Rm oc B had changed. fl2l|. The crossover mag- 
netic field (B* ) between Rm oc B 2 and Rm oc B increased 
with an increase in x, while its absolute value of dRM/dB 
decreased with an increase in x [26j]. 

In order to clarify the contribution of Dirac fermions to 
the electronic transport, we analyzed the results of mag- 
netoresistance and Hall coefficient using a two-carricr- 
type semiclassical approximation. In the low B limit, the 
zero-field resistivities (p(0), Rm and p yx ) are described 



P(P) 
Rm 
p yx 



e(n e p e + n h ph) 

P( B ) ~ /°(°) = n e n h p e p h (p e + p h ) 2 B 2 

p(0) {n e p e + n h p h ) 2 

(-n e pl + n h p 2 h )B 
e(n e p e + n-hPh) 2 



(1) 
(2) 
(3) 



where n c and are the carrier numbers, and p e and 
Ph are the mobilities of electrons and holes, respectively. 
The resulting fact that Rm and p yx are quadratic and 
linear against B in the low B regime is consistent with 
the two-carrier-type semiclassical approximation under a 
low B strength. 

For estimating transport paramatcrs using the above 
three equations, an additional constraint is necessary. 
It has been reported by angle resolved photo emission 
spectroscopy (ARPES) that the Luttinger volumes of 
both electron and hole Fermi surfaces are comparable 
with each other in Ba(FcAs)2 and that Ru sub- 

stitution will not change the Fermi surface significantly 
[30| . Therefore, one can approximate the system as a 
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FIG. 2: (color online) Color contour maps of mobility for electron (p c ) (a) and hole (j/h) (b) in x - T phase diagram estimated 
from two-carrier-type semiclassical approximation. (c) x dependence of p c , p^, and the Hall mobility (pnaii) at 2 K. 



semimetal, that is, n e — is applicable. We systemat- 
ically solved p(0), i?M and p yx by using equations (1) - 
(3) under the condition of n e = nh, where Rm and p yx 
are quadratic and linear against B, respectively. 

Figures 2 (a) and (b) show the color contour maps of 
p e and /ih in the T - x phase diagram. Both p e and ph 
increased with a decrease in temperature in the whole x 
region below the structure and magnetic transition tem- 
perature (T s ) pp| . The enhancement in p c below T s was 
larger than that of p^. At low temperatures, p c and p^ 
decreased as a function of x. The x dependence of mo- 
bility estimated from i?n (mh = Rna.\\/ p(0)) is shown in 
Fig. 2 (c) together with p e and /ij, at 2 K. The value 
of //Hail was significantly smaller than p c and p^- In 
a multi carrier system, the contribution of electron and 
holes on p yx compensates each other. Therefore, the sig- 
nificantly smaller values of /iHaii than those of p c and 
/ih may be indicative of the fact that electrons and holes 
contribute comparably to the electrical conductivity in 
Ba(Fei_ 2 ,Ru a; As)2. 

It appears that the enhancement in mobility below T s 
is related to the occurrence of the Dirac cone, because 
the Dirac cone states form under the SDW order states. 
In the case of the Dirac cone type linear electronic dis- 
persion, the cyclotron effective mass (m*), being propor- 
tional to the inverse of mobility is described in 
Eq. 4 by Fermi energy (Ep), Fermi velocity (vp) and the 
carrier numbers of Dirac cone (jimr) [H- 

to* = Ep/vp = (7™ M r) 2 /wf (4) 
B* = (l/2ehv F ){E F + k B T) 2 (5) 

The experimental observation of 1/p oc bmr 1 ' 2 , 
therefore, indicates low-dissipative transport of Dirac 
fermions, independent of impurity scattering. In a multi- 
band system, the electrical transport is represented by 
mobilities and carrier numbers of all pockets (See Ap- 
pendix). If we approximate both magnctorcsistance and 



the Hall coefficient in the n-carrier system using the two- 
carricr-typc semiclassical model, transport parameters in 
each pocket in the n-carrier system can be merged into 
n e , nh, p e , and ph in the two carrier model. Therefore, 
we cannot distinguish electrons and holes in the two car- 
rier model. Here we thus employed the averaged mobility 
(Mavc = {Pc + Mh)/2) to evaluate the dependence of mo- 
bility on the carrier number of the Dirac cones. In order 
to estimate the carrier number of the Dirac cones, we 
employed a theoretical model for the analysis of magne- 
toresistance in the quantum limit of a Dirac cone [3l| . 
because the carrier number of the Dirac cones cannot be 
estimated directly in the Ba(FeAs)2 system by employ- 
ing the two-carrier-type semiclassical approximation. For 
Ba(Fei_ 2 ,Ru 2: As)2, the evolution of linear magnetoresis- 
tance as a function of B was observed under the SDW 
ordering state. This is consistent with the quantum limit 
of the Dirac cone states [2(| . Considering both the split- 
ting of Landau levels and the thermal excitation energy, 
the onset B (B*) of the quantum limit of a Dirac cone 
can be described by Eq. 5 [3l| . Employing the estimated 
E F and vp from Eq. 5 as shown in Fig. 3 (a) (26|, we 
can estimate timr and to* using Eq. 4. 

For Ba(FeAs)2, it has been found from the first princi- 
ples band calculations that four pockets including two 
Dirac cones with different carrier numbers take place 
in the SDW order state Experimentally, the 

Shubnikov-dc Haas oscillation (ShdH) has shown three 
pockets, however the smallest Dirac cone cannot be ob- 
served [3, while ARPES 0], the de Haas- van Alphen 
effect [8|| and the complex term in THz optical conduc- 
tivity [13j indicate the existence of another small Dirac 
cone. Since the Ep estimated from the magnetoresistancc 
is close to that of the smallest Dirac cone for Ba(FeAs)2 
@i EH j "-mr can reasonably be assigned to the smallest 
Dirac cone. 

Figure 3 (b) shows the dependence of l//i aV c on 



4 




01 

0.05 0.10 0.15 o 0.05 0.10 0.15 

x ( Ru > x (Ru) 



(C) 



35 2 




1 2 3 4 5 6 
„i/2 (1Q5 i/ cm ) 



0.05 0.10 0.15 

x(Ru) 



FIG. 3: (color online) (a) a; dependence of Fermi en- 
ergy (Ef) and Fermi Velocity (vf) of Dirac cone [2(| for 
Ba(Fei_a;Rua;As)2. (b) Carrier numbers of Dirac cone (umr) 
dependence of inverse of averaged mobility (/i ave = (fie + 
A*h)/2) for Ba(Fei_ a: Ru a ;As)2. (c), (d) x dependence of re- 
laxiation time and mean free pass for Ba(Fei_ I Ru a; As)2 at 
2 K. A solid line shows the calculated distance between Ru 
impurities. 



Hmr 1 '' 2 . The value of l/fi avc developed linearly against 
riMR 1 '' 2 . Given the continuous increase in carrier number 
for the two Dirac cones with different carrier numbers, 
both values of m* linearly evolve independently with each 
other as a function of square root of the carrier number of 
the individual Dirac cones. Therefore, even if both Dirac 
cones have a significant contribution on fx ave , l//x ave will 
show a linear evolution as a function of jimr 1 ■ Conse- 
quently, the present observations of l//x ave oc ijmr 1 ' 2 are 
consistent with the expected increase in carrier numbers 
with cyclotron effective mass in the Dirac cone. 

In order to clarify the effect of impurity scattering on 
the Dirac cones, estimation of the mean free path (?) is 
important. Employing both E-p and Dp, which are as- 
sociated with the smallest Dirac cone as shown in Fig. 
3 (a), m* can be calculated from Eq. 4. We accord- 
ingly estimated both the relaxation time (t) and I using 
Mavo, m* and vf- Note that both t and I for the smallest 
Dirac cone may be underestimated because ^ avo is also 
influenced by other larger pockets. Figure 3 (c) and (d) 
show the dependencies of both r and I as a function of 
x. Both r and I are almost constant with x, and I is 
larger than the average distance between Ru impurities 
above x = 0.05. Since I is underestimated for the small- 
est Dirac cone, the present estimation indicates that I 
is not reduced by Ru substitution. Therefore, the ob- 
served almost constant value of I, independent of the av- 



erage distance between Ru impurities, is consistent with 
the expected large suppression of backward scattering of 
Dirac fcrmions. 



IV. SUMMARY 

We studied electronic transport of Dirac fermions in 
iron pnictide materials Ba^ei^Rur As)2 using rapidly 
quenched single crystals. The electronic transport of 
parabolic bands was suppressed by scattering due to the 
crystal lattice distortion as well as the Ru impurity, while 
that of the Dirac cone is not significantly influenced. 
From the analyses of magnetoresistance and Hall coef- 
ficient using the two-carrier-type semiclassical approxi- 
mation, it was found that the inverse of the average mo- 
bility, being proportional to the cyclotron effective mass, 
developed as square root of carrier number of the smallest 
Dirac cone. This can be understood in terms of the in- 
trinsic low-dissipative carrier transport of Dirac fermions. 
The estimated mean free path, using the effective mass 
of the smallest Dirac cone [2(| and the average mobility 
obtained from the semiclassical analysis, was constant for 
x and larger than the average distance between Ru impu- 
rities. This is consistent with large suppression of back- 
ward scattering in Dirac fermions. The present results 
have clearly demonstrated that the low dissipative carrier 
transport as a consequence of the nature of the berry's 
phase of the Dirac cone is realized in Ba(Fei_ a; Ru a ;As)2 
under the SDW ordering state. This suggests that a pos- 
sible enhancement of the phase coherence of Cooper pairs 
via the Dirac cone high mobility state plays an important 
role for missing in temperature fluctuations in the pnic- 
tide family. 
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Appendix: Magnetoresistance and Hall coefficient 
for 4-carrier-type semiclassical approximation 

It was theoretically predicted by the band calcula- 
tions employing the local density approximation and 
the density functional theory that the Fermi surface of 
Ba(FcAs) 2 is composed of 4 pockets including 2 Dirac 
cones with different carrier numbers in the SDW order- 
ing state 0, HI]. In this case, both magnetoresistance 



and Hall coefficients can be described using mobilities 
and carrier numbers of four pockets. Here, we showed 
equations of magnctorcsistance (Rm) and Hall resistiv- 
ity (p yx ) derived from the 4-carrier-type semiclassical ap- 
proximation in the low-B limit. 

A semiclassical-type conductivity tensor of a carrier 
can be described in eq. (A.l) (33j . 



en i p l 



i + ni 



1 ihB 

-ihB i 



(A.l) 



In the 4-carrier-type expression, the resistivity tensor 
can be obtained using the inverse of sum of 4 conductivity 
tensors. 



Pxx Pyx 
Pyx Pxx 



(A.2) 



In the low-£> limit both Rm and p yx can be approxi- 
mated as shown in eqs. (A. 3) and (A. 4). 



Rm 

Pyx 
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P = (( V MiMl^i - V MiM3"i) 



7 
<5 
c 



+(\/ M2M|"2 - \J plp3,n 2 ) 2 )ns, 
{\J L L il4. n i ~ \j p\P2n\) 2 n 2 

(Hini + P2n 2 + P3m + M4"4) 2 

Pyri\ + p\n 2 + nl~n 3 + p\ni. 



(A.3) 
(A.4) 



where 

= ((\/ 1 PiP\n\ - \J pIpa^) 2 + (y ' p 2 p\n 2 - \j plp A n 2 ) 
+(\/M3M4 Tl 3 - \J plPin 5 ) 2 )n A 
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